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interfering paths 



M. Ozana and A. Shelankov* 



o 
o 

(N 

Or 
<D 

in 

(N 



S3 
O 
O 

Oh' 
o3 



13 
G 
O 

o 



> 

oo 

On 
O 



c3 



i 

C 

o 
o 



X 



Department of Theoretical Physics, limed University, 901 87 Umed, Sweden 
(Dated: February 1, 2008) 

We apply the method of two-point quasiclassical Green's function to geometries where the trajec- 
tories include interfering paths and loops. For a system of two superconducting layers separated by 
partially transparent interface, corrections to the quasiclassical solutions for the Green's function 
are explicitly found as well as the deviation from the normalization condition. 
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The method of quasiclassical Green's function has been 
long known as the most efficient tool in the theory of 
superconductivity. In the quasiclassical theory, the su- 
perconductor is described by the matrix Green's func- 
tion which depends on the coordinate and the momen- 
tum as well as .the energy variable. It obeys the 
berger equation^ or its nonstationary generalisation! 
the form of which is similar to the Boltzmann kinetic 
equation, with the distribution function replaced by a 
matrix. Following the analogy between the Eilenberger 
and Boltzmann equations, one may say that the Eilen- 
berger equation describes particles moving along classical 
trajectories, the matrix nature of the "distribution func- 
tion" accounting for the internal quantum electron-hole 
degree of freedom of the quasiparticle in a superconduc- 
tor. As in the Boltzmann equation case, the notion of a 
smooth trajectory remains meaningful in the presence of 
disorder: Although individual impurities may lead to a 
large angle scattering, the averaging with respect to the 
impurity positions restores the spatial homogeneity, and 
the trajectories remain smooth. The presence of disor- 
der (leaving the localization effects out of the picture) is 
seen only in a finite time that the particle lives on a given 
trajectory. 

In a mesoscopic superconductor, where during its life 
time the particle may suffer multipe reflections on the 
boundary, or a bulk sample with interfaces, the typical 
trajectory may be rather complicated. Compared with 
classical billiards, trajectories in a superconductor be- 
come even more intricate due to the Ajadreev reflection 
processes. It has been recently argued inQ that the deriva- 
tion of the quasiclassical equation is heavily based on the 
assumption that the trajectories are single-connected i.e., 
there exists the unique path which connects any pair of 
points visited by the particle: The one dimensional char- 
acter of the trajectory allows one to single out the phase 
factor due to the motion with the Fermi momentum and 
derive an equation for a slowly varying envelope, which 
is the Eilenberger equation. Otherwise, if the particle 
has alternative paths to go from one point to another, 
the quantum interference between the paths takes place, 
which is sensitive to variations of the path on the scale 
of the Fermi wave length Xp. In this case, the Fermi 



wave length oscillations cannot be disentangled from the 
behaviour at the larger spatial scales of order of the co- 
herence length. Although interfering paths always exist 
in any realistic situation (e.g. due multiple scattering 
on impurity), the interference contribution is usually in- 
significant since it vanishes after averaging with respect 
to the parameters of the trajectories. This is the reason 
why the quasiclassical Eilenberger or Boltzmann equa- 
tions, which completely ignore the interference, work very 
well. 

However, there may be special reasons for the interfer- 
ence to survive the averaging. For instance, the interfer- 
ence of waves propagating along two time reversed path 
is always constructive if time reversal is a valid symme- 
try, and it leads to the well-known phenomena of weak 
localization. The interference may be robust in certain 
billiard-like geometries where the alternative paths are 
created by reflections on the boundaries or interfaces. 
Here, one may invoke the analogy with a classical billiard 
where intergrability is very sensitive to the geometry. As 
an example, consider an ideal two-layer system shown in 
Fig. |l|. Trajectories for the case of equal thicknesses of 
the layers is shown schematically in Fig. |^ (a). For any 
direction of incidence, the two out-going waves created 
by the collision with the specular interface combine to- 
gether after reflection by the outer walls. The trajectory, 
that is, the set of the points spanned by the particle, is 
not single-connected. The lengths of the two alternative 
paths are the same so that their interference is insensi- 
tive to the angular averaging and, therefore, may be of 
importance. (Whether the interference indeed influences 
observables requires an additional analysis.) The simple 
of the two recombining paths picture Fig. ^| (a) exist 
only if iwht a high precision the layers have equal thick- 
nesses, a® = a' r ). However, the interference is possible 
for any relation between aW) and a^: For the arbitrary 
parameter a^/a^"\ rational or irrational, the wave frag- 
ments (partially) recombine after two collisions with the 
outer walls (see Fig. §(b)). It is only when the assump- 
tion about the specular character of reflection is relaxed, 
it becomes statistically impossible for the partial waves 
created by a collision with the interface to meet each 
other ever again. In these conditions when the alter- 
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native paths are absent, the trajectory becomes simply 
connected, or tree-like in the terminology offl Using the 
method developed inH, one is able to calculate the quasi- 
classical Green's function on the tree-like trajectory. 
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FIG. 1: Sandwich build of two superconductor layers with 
the order parameters A' !,r ' and thicknesses a' i,r ' separated 
by a partially transparent interface. 
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FIG. 2: Traejectories in a double-layer system with a par- 
tially transparent interface. In the ideal system, the reflec- 
tions are specular, (a) When thicknesses of the layers are 
equal, the trajectories linked by the collsion with the inter- 
face at the point A meet each other again at the point B after 
the specular reflections by the boundaries, (b) For a general 
irrational relation between the layer thicknesses, the trajec- 
tories meet each other meet again after two reflections. The 
lengths of the two interfering paths connecting the points A 
and B are equal for any direction of the trajectory coming to 
the point A so that the interference survives the momentum 
averaging. 

The putpose of this paper is (i) to substantiate the 
conjectured about the qualitative importance of the 
single-connectivity of the trajectories in the conventional 
quasiclassical theory; (ii) to suggest a method which al- 
lows one to include the possibility of interfering paths 
in the quasiclassical theory of superconductivity. Being 
interested mainly in general questions of the quasiclassi- 
cal Green's function theory, we consider for definiteness 
the simplest possible system where the interfering path 
exist, namely, a sandwich shown in Fig.|l|. For the spec- 
ular case, this system is simple enough to solve the full 
Gor'kov equation and thus to learn "exact" properties. 
The exact solution being known, one can judge on the va- 
lidity of approximate schemes and analyse the physical 
reason for their failure-df it happens. It has been demon- 
strated numerically incl that the exact and quasiclassical 
theory does give a different result in this case. To gain a 




FIG. 3: Trajectories in a sandwich is half-transparent inter- 
face. The intensity is proportional to the logarithm of the 
absolute value of the exact normal metal Green function after 
Gaussian average in a small region of the momentum parallel 
to the interface. 



better qualitative understanding of the numerical results 
ofl is the main goal of this study. 

It is important to realize that the exact solution to the 
Gor'kov equation must not be directly compared with the 
quasiclassical solution. For instance, it is clear that the 
momentum resolved density of states (DOS) found from 
the full Gor'kov equation for a layer with specular walls, 
shows (5-function peaks as a function of energy or mo- 
mentum due to the space quantisation in the layer. On 
the other hand, the quasiclassical DOS is smooth as in 
an infinite layer. As discussed irfl, it is a coarse-grained 
Gor'kov Green's function that should be compared with 
the quasiclassical results. The coarse-grained function is 
generated by integration of the full Green's function in 
a small region of x — p space. In particular, one regains 
a classical picture of trajectories on the coarse-grained 
level of description. As an illustration of this point, 
Fig. U shows the picture of trajectories produced by the 
integration of the exact normal metal Green's function 
G(x, exp[ip||(r — r')] in a small region of pn, p|| 

being the in-plane momentum. 

The quasiclassical theory of superconductors with in- 
terfaces and walls has been discussed in many papers. 
Surfaces and interfaces violate the condition of applica- 
bility of the quasiclassical theory, and their presence is in- 
corporated in the quasiclassical theory through a bound- 
ary condition, again similar to the Boltzmann trans- 
port equation. For the case of an impenetrable wall, 
the boundary condition has been derived in Ref.Q. The 
matching condition for a partially .transparent suecular 
interface has been derived in RefH as well asErEl. An 
isolated interface does not pose severe problems for the 
quasiclassical theory of superconductivity. The case of 
two and more interfaces turns out to be more difficult. 
Analysing the solution to the full Gor'kov equation for 
a sandwich made of two superconductors with a pas- 
tially transparent interface, some authors (seeE3, andtll) 
have come to the conclusion that the normalization of 
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the Green's function is violated and, therefore, the qua- 
siclassical scheme fails in this case. To a large extent, 
the present paper has been motivated by the observation 
made inEJ. 

The plan of the paper is the following. In SectionJL we 
review the method of the 2-point Green's functionljQ on 
an infinite trajectory. In Section [n| we derive quasiclassi- 
cal Green's function for an isolated loop, and discuss the 
normalisation condition. We show that the normalization 
condition is violated before the coarse averaging, but not 
after. In Section [II, we consider two coupled loops, and 



energy, 



show that the normalisation condition is not valid even on 
the coarse-grained level. Finally, in Section [V, we apply 



the results to a superconducting two-layer and analyse 
applicability of the quasiclassical theory. The physical 
intepretation of the results is given in Section^. 



I. 2-POINT QUASICLASSICAL GREEN'S 
FUNCTIONS 

We start with a short review of the formulation of the 
quasiclassical theory in terms oiuthe two-point Green's 
function on classical trajectoriesaO. Here, we present a 
novel derivation of the equations in a form convenient for 
further generalisations. 

The main assumption of the quasiclassical theory of 
superconductivity is the existence of two distinct spatial 
scales: the short one is the Fermi wave length Xf = 1/??f, 
Pf being the Fermi momentum, whereas the supercon- 
ducting coherence length, £o, exemplifies the large scale. 
The quasiclassical theory is valid in the zero appjopdma- 
tion with respect to the small parameter Ai?/£o™^- In 
effect, the Fermi wave length is put to zero and the orbital 
motion is essentially classical, that is, the particle coor- 
dinate belongs to a chosen classical trajectory. In this 
picture, the electron-hole mixing, which is the essence of 
superconductivity, is a quantum internal degree of free- 
dom. 

For a given momentum, the direction of velocity is op- 
posite for the electron and hole components of the wave 
functions, and the derivation should start in a way that 
both direction on the trajectory enter on equal footing. 
Consider propagation on a straight line and choose the 
x— axis along the line. Denote the matrix 2x2 propaga- 
tor &(x\,X2)- Throughout the paper, we deal with the 
retarded Green's function of the Keldysh technique, G R ; 
to simplify notation we omit the superscript R. The sta- 
tionary Green's function satisfies the Gor'kov equation 



(e-W)G(x 1 ,x 2 ) = ltf(xi 



X 2 ) 



(1.1) 



where e = e + iS is the energy variable of the retarded 
Green's function, and the matrix operator H reads 



£ A 
A* -£ 



(1.2) 



where A(x) is the superconducting order parameter, S 
is the matrix (retarded) self-energy, and £ is the kinetic 



Px 



Pi 



2m 2m 



Px 



1 d 

i dx ' 



with the Fermi energy subtracted. For simplicity, we 
consider the isotropic case in the absence of a magnetic 
field. It is straightforward to take into consideration the 
anisotropy of the order parameter and the vector poten- 
tial. 



The boundary condition to Eq. (1.1) is 

G(xi,x 2 ) — * , |xi-x 2 |^oo 



(1.3) 



as it follows from the requirement of the analyticity of 
the Green's function. 

Considered as a functio n of the first argument x\, the 
Green's function G in Eq. (1.1) describes the wave propa- 
gating from the point x 2 ■ If the energy variable e is small 
in comparison with the Fermi energy p F /2m 1 the wave 
has two components with the momenta close to +pf or 
— pf- In the quasiclassical theory, these components are 
ascribed to two trajectories with the opposite momen- 
tum, and the corresponding trajectory propagators obey 
separate inhomogeneous equations. To carry out the sep- 
aration, we use the following procedure. 



The 5— function inhomogeneity in Eq. (1.1), that is the 
source of the wave, is equivalent to the following condi- 
tions at the points x± = x 2 ± 



(xi,x 2 ) 



| X1—X2+O 

X\—X2~0 
£1=212+0 







1 d ^, 

G(xi,x 2 j| 
2mdx\ lxi=x 2 -Q 

Defining the differential operators D^, 



(1.4) 
(1.5) 



D , 



1 1 

± 



2 2ipp dx 



one rewrites Eqs. ( |l.4| ) and (1.5), as the following two 
inhomogeneous equations, 

WfD G{Xi,X2)\ = 1 . 

Ia?l= X2— 

These conditions can be satisfied by the following de- 
composition: 



G(xi,x 2 ) = G <T (xi,x 2 ) 



(1.6) 



where the partial pro pag ators obey the homogeneous 
linear equation Eq. (1.1) at x± ^ x% and the inhomoge- 
neous conditions 

iv F Dl i G' J {x 1 ,x 2 )T z \ x ^ + ' J - a =<Wi , a,a' = ±. 

\xi—X2—cr-0 

(1.7) 

which replace the (S-function source in the r.h.s. of 



Eq. (1.1) 
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Now, one observes that <G + and G~ correspond to the 
opposite direction of the Fermi momentum i. e. they can 
be presented as 



z (x 1 ,x 2 )t z = g ± (x u x 2 ) e =MM*i-«*) , (1. 



where g^ (xi , x 2 ) is a slowly varying function if the energy 
e is small in comparison with the Fermi energy; the factor 
with the Fermi velocity vf = Pf /tti and the Pauli matrix 
t z are introduced for future convenience. 
Indeed, for slowly varying g^, 

-4-g ± {xi,x 2 )e ±lpF{xi - X2) « ±p F g ± (x l ,x 2 )e ±lpF(xi - X2) 
i dx\ 

(1.9) 

Assuming that typically g ± {xi 1 x 2 ) vary on the coherence 
length £oj the corrections to this relation are of order of 
^f/£o- Small quantities of this order are consistently 
neglected in the quasiclassical theory. Note that-this for- 
mula correspond to the Andreev approximations where 
the derivatives of quasiclassical envelope functions are al- 
ways neglected if they stand next to pp. 



With the accuracy Eq. fll.9|) holds, are projectors, 



D" 



<t'=± 



(x l7 x 2 ) = G a (xi,x 2 ) 



The projecting pro perty guarantees that one of the con- 
ditions in Eq. (|L7|), namely, corresponding to a' = —a, 
is automatically satisfied. The remaining conditions read 



g + (x ll x 2 )\ 



g (351,2:2) 



I xi— X2+0 
I X1—X2—O 
I xi— X2— 
I x 1 — x 2 +0 



(1.10) 

(1.11) 



It is clear that g + and g~ are associated with the Fermi 
surface points +ppn and —ppn, respectively, n being the 
orientation of the x-axis. Spanning different orientation 
.of the x— axis, one defines function g[x\,x 2 \n,R) where 
n and R are the parameters of the trajectory, and X\ i2 
are the coordinates on the trajectory see Fig.|. Then in- 
formation on the direction of the momentum is contained 
in n, and the superscript ± in g , being redundant, can 
be omitted. 



Plugging the ansatz Eq. (1.8) into Eq. (1.1), and using 
the Andreev approximation 



-IVp 



d_ 

dx 



one obtains the first order differential equation for g. 

The 2-point Green's function g e {x\, x 2 \n, R) obeys the 
following set of equations^: 



iv F-Q^ $(71) +H ei „.(ri) ) g s (xi,x 2 \n,R) = iv F 5(x 1 -x 2 ), 



: (x 1 ,x 2 \n,R) ( -iv F — <J>(r 2 ) + H E , n (r 2 ) ) = iv F 5{xi - x 2 ), 



(1.12a) 
(1.12b) 



where the derivative in Eq. ( 1.12b ) operates backwards. 
The (5-functions i n the r.h.s. generates the discontinuity 
required by Eqs. 



1.10) 



In Eq. (1.12), $ is the gauge invariant scalar potential, 



and H E n is 2 x 2 traceless matrix, 



E-Vp p s 

A* 



-An 

-e + vp ■ p s 



Vp 



where A n is the order parameter (which may depend on 
the direction n), and p s = —-A, A being the vector 

potential, and S is built of the impurity self-energy as 
well as the part of the electron-phonon self-energy not 
included to the self-consistent field A. Below, we omit 
R, n and e for brevity. 



The boundary co nditi on to Eq. (1.12) on an infinite 
trajectory from Eq. (|l. 



g(xi,x 2 )—>0 for \xi-x 2 \ 



00 



(1.13) 



The equations of motion Eq. |4yj) can be solved by 
factorizing the Green's functionEIO. For this, one finds 
linear independent two-component solutions 4>±(x) = 
(« ± ) *° the equation 



iVF-5- + H 
ox 



= 



4>±(x) — > , x — > ±00 
(1.14) 



where H is in the operator on the l.h.s. of Eqs. (1.12); 
This equation is essentially the same as the one first de- 
rived by Andreev in the theory of reflection on the NS i in i 
terface. The ratio u±/v± obeys the Riccati equationtilEJ, 
which is most convenient for numerics. Known solution 
to the Riccati equation for u±/v±, one readily finds full 
<f>± as explained inO. 



The solutions to the linear equation Eq. (1.14) can al- 
ways be normalized as 



(1.15) 
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FIG. 4: Classical trajectory. When no interfaces or bound- 
aries are present the trajectory is a straight line. It can be 
specified by an initial point R and the unit vector n, which 
shows the direction of the trajectory. Introducing the coor- 
dinate along the trajectory x, any point r on the trajectory 
can be written as r — R + xn. 



where here and below the bar on the top of a 2-component 
column, tp, or a matrix X, has the following meaning 



1 



V> =tP Ty- 



X 



TyX Ty 



(1.16) 



where T stands for transposition and r y is the Pauli ma- 
trix. Taking advantage of the property H = —H, one 
can prove that the norma lizati on in Eq. (1.15) is gener- 
ally compatible with Eq. (1.14). 
Denote also 



Then, 



4 = H-v 



<>v$v> = <W 



± 



^2 $ v< & 



(1.17) 



as immediately follows fro m the normalization Eq. ( 1.15| ) 
and the definitions in Eq. ( 1.16 V 

Now, the solution to Eq. (1.12) reads 



g{xi,x 2 ) = 



<j)+{xi) 4> X + {x 2 ) 

(j)-{xi) <j) t _{X2) 



xi > x 2 ; 
Xi < x 2 ■ 



(1.18) 



In the standard formulation of the quasiclassical tech- 
nique after Eilenberger, one deals with the 1-point 
Green's function Site), which is related to the 2-point 
Green's function as£3 



.(*) 



(1.19) 



where 



9±{x) = g(xi =x±0, x 2 = xtO) 



Combining Eqs. (1.12a), and ( 1.1,2b ), one re-derives the 
well-known Eilenberger equationtU, 

iv F • Vg + Hg -gH = 0, 



for the matrix function g. The expressions for observ- 
ables like the charge current and the density of states are 
presented in Section |X| 



In the factorized representation Eq. (1.18), the normal- 
ization condition for the Eilenberger function Eq. ( 1.19 ) 



(1.20) 



imm ediat ely follows from the orthonormality relations in 
Eq. (p7|) . 

To trac e the general features needed for the validity of 
Eq. ( 1.20| ), we first note that 



9+{x) - g-(x) = 1 



as is dictated by the 5— function in the r.h. s. of Eq. (1.12). 
Taking square of this equation a nd E q. ( 1.19 ), one gets 
the normalization condition Eq. (1.20) provided 



g + (x)g-(x) = g-{x)g + (x) = . 
The proof of these relations present ed in Ref.El, essen- 



tially uses the asymptotics in Eq. ( 1.13 ), which makes 
sense only for an infinite trajectory. Besides, the proof 
tacitly presumes that the quasiclassical equations of mo- 



tion Eq. (1.12) are applicable everywhere on the trajec- 
tory. The latter assumption is not valid if the trajec- 
tory under consideration is crossed by another trajec- 
tory - the trajectory crossing happens, for instance, on 
partially transparent interfaces. As shown in Reffl, the 
trajectory crossings do not influence the normalization 
provided the crossings do not change the simple one- 
dimensional topology of the trajectories. It is the purpose 
of the present paper to analyse more general cases. 

After the review of the two-point quasiclassical Green's 
function method, we will extend the method to the case 
when the classical trajectories have more complicated 
topology. In the standard theory, the Fermi wave length 
scale doe s not enter the theory after the factorization in 
Eq. (1.8) has been carried out. As we will see later, this 
disentanglement of the different scales can not be per- 
formed in the case of more complicated geometry. 



II. GREEN'S FUNCTION ON A LOOP 

The theory reviewed in the previous section refers to 
the situation when the trajectory is a straight line. One 
can easily generalise the derivation to the case when the 
particle moves in a smooth quasiclassical potential and its 
trajectories are not straight lines. If the trajectory cor- 
responds to the infinite motion, no impo rtant changes in 
the theory occur: one can use Eq. (1.12) with the under- 
standing that x is the coordinate along the path. If the 
orbit is closed, the situation is very different. Physically, 
the obvious difference is that if the particle moves along a 
trajectory it eventually comes again to the starting point; 
then, statements like x\ > x 2 do not have unequivocal 
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sense. Also, the boundary condition in Eg . (1.13), with- 
out which the equation of motion Eq. ( 1.12 ) do not have 
a unique solution, becomes meaningless. 

Of course, in a clean homogeneous metal one can 
hardly imagine any smooth loop-like trajectory because 
any potential is screened by mobile charges. Actually, 
we have in mind to consider loops created as a result of 
multiple reflections by interfaces. However, for the sake 
of presentation, it is convenient to us to analyse first a 
smooth closed trajectory. Besides, the problem is of gen- 
eral interest. 

Consider a closed loop trajectory sho wn in Fig||. We 
accept it without derivation that Eq. (1.12) may serve 
as the equation of motion for the two point trajectory 
Green's function. What is needed is the boundary con- 
dition replacing Eq. ( 1.13| ). 

First, we recall that the full propagator is given by G 
in Eq. (O), 



iv F G(x u x 2 )t z = g(x u x 2 ) e **( ai - aa > , 



(2.1) 



omitting the redundant superscripts ±. Consider 
the Green's function at nearly coinciding arguments, 
G(x± = X2 + 0,2:2). Using equation of motion, one can 
find G(x±, X2) advancing x\ in the direction of arrows in 
Fig.||. Eventually, one gets again close to X2 when x\ 
reaches the value x 2 + C — 0, C being the circumference 
of the loop. On the other hand, the point x\ = X2 + £ — 
is the same physical point as x\ — X2 — 0. From the re- 
quirement that the Green's function is uniquely defined, 
one gets the following boundary condition 

G(xi = x 2 + C, x 2 ) = G(xi = x 2 - 0, x 2 ) ■ 

From here one gets the following boundary condition 
for the quasiclassical Green's function 



z ipFC g{x 2 + C, x 2 ) = g(x 2 - 0, x 2 ) 



(2.2) 



where g(x2 + C,x 2 ) is the quasic lassic al Green's function 
continued with the help of Eq. (1.12) from x\ = X2 + 
along the loop of the length L to the point x\ = x 2 + C. 

X] x 2 



As the next step, we derive formally exact expres- 
sion for the 1-point quasiclassical Green's function via 
the quasiclassical evolution operator U(x. x'). The latter 
solves the homogeneous part of Eq. (1.12), 



■1 + h) 



ivftt- + ti I U(x, x') = 
ox J 



(2.3) 



with the initial condition U{x, x) = 1. A formal solution 
to Eq. (2.3) reads 



U(x, x 1 ) = exp 



1 



dx H 



(2.4) 



where x— ordering of the matrices H(x) is implied. Pr< 
erties of U has been discussed in many papers, see e.g. 111 
In particular, det {7 = 1, and 



U(x,x')U(x,x') = 1 



(2.5) 



The expression for the evolution operator in the Riccati 
equation technique is presented in Section [b] (see alsoQ) . 
Denote 



U x = U(x + C, x) 



U T . = e ippC U x 



Then g(x + C, x ) — U x g(x + 0, x) and the boundary con- 
dition Eq. (2.2) acquires the form 



U x g{x + 0, x) = g(x - 0, x) 

This equation together with 

g(x + 0, x) — g(x — 0, x) = 1 

allows one to obtain the expression for the one-point 
Green's functions via the evolution operator: 



9+(x) = 
9-i x ) = 



1-Z4 

u x 
\-u x 

i-u x 



(2.6a) 



(2.6b) 



For illustration, we consider some simple applications of 
this result. 



Homogeneous loop 



FIG. 5: A loop. The Green's function is a function of the two 
coordinates along the loop Xi and X2- As shown in the picture, 
advancing gradually the first argument from the position Xi = 
X2 + O in the direction of its increase, one approaches the point 
marked by the dashed circle which locally correponds to the 
coordinate x\ = X2 — 0. 



Suppose the "Hamiltonian" H in Eq. (2.3) does not 
depend on the coordinate. T he evolution operator U x 
can be easily found from Eq. (2.4). For this, we present 
H as 

H = £ Ml ~ £ V-Vi 
where ip± are the eigenfunction of H, 

ffV±=±^± , »£>0 (2.7) 
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nor maliz ed in accordance with Eq. (1.17) 
Eq. (2.4), the loop evolution operator reads 



From <7 + Eq. (2.6a), 



u x =u tjj+^ l + + - i>-i>\ 

u 

independent on x. Parameter u, 



u = e v F 



(2.8) 



(2.9) 



shows the variation of the phase of the wave and its abso- 
lute value upon propagation around the loop. The mod- 
ulus of u, \u\ = exp[— L Ss^/vf] < 1, is less t han unity. 
The one-point Green's function Eq. ( |2.6b| ) reads 



1 
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u e 



ip F C 



-ipF C 



1 -ue lpFl 



ip+w+ - - — 
i — 



-ip F C 



ip-^t . (2.10) 



This result is rather transparent physically. The 
Green's function has a contribution from (quasi)electrons 
and (quasi) holes given by the first and second terms, re- 
spectively. The (retarded) Green's function is built of 
waves decaying in the direction of propagation, and one 
can use attenuation of the wave, as a meter which 
shows the length of the path the particle has travelled. 
Expanding the Green's function in Eq. ( 2.10| ) with re- 
spect to u, 

g = ijj + vb\ (l + ue lpFC + («e' Pf£ ) 2 + . 

- V-V'l (i + ue- iPFC + («r' Pf£ ) 2 + ...) . 

one gets information on the paths contributing the prop- 
agator and their length. Looking at the first line, one 
sees that the Green's function has a contribution which 
corresponds to the electron-like quasiparticle ip + circling 
the loop in the positive direction (where the phase in 
Fermi momentum factor increases); the paths have the 
winding numbers 0,1,2... and the lengths 0, C, 2C, . . . , 
respectively. The second line shows the contribution due 
to holes: Decaying in the direction of propagation, they 
circle the loop in the negative directipjs as can be read 
off the Fermi momentum phase factored. In the limit of 
an infinite loop, when u — > 0, one restores the well-know 
result for a homogeneous superconductor. 

The electron and hole parts in Eq. ( 2.10] ) have poles 
at the energies where the denominators have zeroth. In 
the ballistic case, when the non-Hermitian part of H is 
only due to the infinitesimal 5 in e + iS, the poles give 
the bound states which correspond to the usual Bohr- 
Sommerfcld quantisation for closed orbits. Impurities, 
generating the imaginary part of the self-energy, reduce 
u and broaden the bound states. For quantities averaged 
with respect to disorder, the closed topology of the orbit 
is seen only if the orbit is shorter than the mean free 
path. At Matsubara frequencies e = iu) n , the orbit must 
be shorter than the coherence length vf/^u- 

A qualitative difference between the cases of a loop and 
an infinite trajectory, is clearly seen in the expansion of 



g+ = ^4( 1 + MeiPF£ + ( Me ' PF£ ) 2 + -) 

- V-V'l (ue- lpFC + (we-^ £ ) 2 + ...) (.2.11) 



One sees that, unlike the case of an infinite trajectory, the 
hole contributes to g + . For this, the hole must circle the 
loop at least once. The reason for the difference between 
electron and hole is physically clear. Indeed, by defini- 
tion, g + (x) is the propagator g(x + 0, x) from the source 
at the point x to the n eighbo uring point x+0. The source 
(i.e. the r.h.s of Eqs. ( 1.12a| ) creates electron which prop- 
agates in the positive direction, and the hole moving in 
the negative direction, nominally, away from the point 
x + 0. However, having made a full turn(s) around the 
loop, the hole reaches the point x + 0, and, therefore, 
does contribute to g + . Ultimately, these "wrong" pro- 
cesses where the excitation effectively moves in the wrong 
direction, lead to the violation of the normalization con- 
dition for the quasiclassical Green's function. 



2. Normalization condition 



To check the normalization condition, we note that 



g — 1 can be conveniently presented as 
9 2 - 1 = 4 g+g- 



From Eq. (2.6), g+g- = U x /(\ — U x ) 2 , or in the eigen- 
function representation, 



aiPFC 



9+9- = 



(1 -ae^) 2 



-ip F C 



(1 - ue^" £ ) 2 



(2.12) 



One can see that the orthogonality and normalization 
condition do not hold for the Green's function on a loop. 
Of course, this is not informative from the point of view 
of checking the validity of the quasiclassical theory which 
deals with co arse- grained objects. Note, however, the 
r.h.s. of Eq. (2.12) vanishes when u — * 0, restoring the 
normalization Eq. (1.20) in the limit of an infinite loop. 



3. Coarse-grained average 

Observables like the local density of states (DOS ) cal - 
culated with the help of the Green's function in Eq. ( 2.6b ) 
or Eq. (2.10), may show wild fluctuations as a function of 
position or energy. Indeed, in a ballistic system the DOS 
is given by a series of closely spaced <5-function peaks due 
to the space quantization in the loop. On a low-spatial 
resolution, when loops with slightly different C contribute 
to observables, DOS becomes an essentially smooth func- 
tion. The summation over the loops in a small interval 
of SC is equivalent to averaging witturespect to the phase 
of tp — pf£ provided pf 5C > 1 (cfllH). 
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In the representation Eq. (2.10) the integration is eas- 
ily performed and the coarse-grained Green's function 

(jcoarse IS simply 



3c 



(2-13) 



For a homogeniuos loop, this expression coincides with 
that for an infinite trajectory. 

This result is readily generalised to the inhomogeneous 
loop. Indeed, the ev olution matrix can be always pre- 
sented as in Eq. (2.8) because it is the only form com- 
patible with Eq. (2.5). Unlike the homogeneous case, ip± 
are cc-dependent and there is no simple way to find them. 
Operationally, ip± are eigenfunctions of U x and u and 1 /u 
its eigenvalues (the latter are global x-independent val- 
ues). Without lost of generality, \u\ < 1 and ip + is the 
eigenfunction corresponding to the s malles t eigenvalue 
of U x . Then, again one comes to Eq. Of course, 

the effective wave functions ip± depend on the coordi- 
nate as well as the length the loop C The normalization 
cond ition is granted by the orthonormality relation in 
Eq. Cl7b . 

We see that the coa rse-gr ained averaging eliminates 
"wrong" paths in Eq. ( 2.11 ) and the averaged Green's 
function can be found from the standard quasiclassical 
theory. However, as we see later, alternative paths may 
survive the averaging in more complicated geometries. 



III. DOUBLE LOOP 

In this section, we generalise the theory to the case of 
two loops tied by a knot, see Fig. ^. 




FIG. 6: Two coupled loops. The arrows show the direction in 
which the Fermi phase increases (the momentum direction). 
The filled circle stands for the knot which ties the loops al- 
lowing for the coupling of the loops. Fhe amplitudes of the 
waves at the points adjacent to the knot, 1, V on the left loop 
and 2, 2' on the righ t lo op, are related by the .S'-matrix in 
accordance with Eq. (3.2). 



As before, we denote the propagator g{x\, x-z) with the 
understanding that x\ and Xi may stand for either or 
x^ r \ the coordinates of the points on the left and right 



loops, respectively. The initial and final points of the 
propagator may be on the same loop or belong to the 
different loops, see Fig. |[ The coordinates x^ 1 '^ are 
counted from the knot in the direction of the correspond- 
ing arrow. They assume values < x^ < C^ l ' r \ C 1 * 1 '^ 
being the circumference of the left or right loop. The 
equa tions of motion within each of the loops are given by 
Eq. ( 1.12| ) . We couple the loops with a phenomenological 
unitary S'-matrix, 



S = 



P t 
t p 



(3.1) 



where parameter t describes the inter-loop coupling 
( "transmission amplitude" ) and p is the intra- loop scat- 
tering amplitude. These parameters are taken same for 
the electron and hole components of the wave function. 
Apart from the the unitarity condition, 



\p\ 2 + W = l 



pt* +p *t = Q, 



p and t are free parameters. In the limit t — 0, one has 
two independent loops, and, if r = 0, the two loops are 
equivalent to a single loop (with the shape like "8" ) . 

Denote ipi,i>, 2,2' the (two-components) wave functions 
at the corresponding point 1,1', 2 or 2', adjacent to the 
knot coupling the loops (see Fig j^) . Then, the S'-matrix 
expresses the amplitudes of the outgoing waves, V'i'.2', 



1p2> 



P t 
t p 



(3.2) 



via the incoming wave amplitude s ip i,2- 

The formal solution in Eq. (2.6) for the one-point 
Green's functions g±{x^ l,rS> ) remains valid if the loop evo- 
lution operator U x (i,r) properly accounts for the knot. 

Consider for definiteness the right loop. The full loop 
evolution matrix can be written as the product, 



) U ( - r \x^ r \Q) K (r) f/ (r) (£ (r) ,a; (r) ) 



(3.3) 



where U^ r \x^ r \ x^ ') is the in-loop evolution matrix 
found from Eq. (2.2), and TZ^ r ' (the transfer matrix) ac- 



counts for the propagation through the knot from the 
point 2 to the point 2' (see Fig . |q ) . The latter can be 
evaluated with the help of Eq. ( |3.2j ) by the following ar- 
guments. 

While calculating the one-point Green's function for 
the right loop, the waves in the left loop obey homoge- 
niuos equation of motion. Then, the amplitudes ipi (see 
Fig.0) and i/jy are related as 



where 



V>1 = U {V) i> v 



u (i) = e i PF cWjj{l) 



(3.4) 



(3.5) 
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and is the quasiclassical evolution matrix in the left 
loop which advances x^ from x"' = to x^ 1 ' = CS- 1 ': 

U® = U«\C®,0). 



Combining Eqs. (|3.2[ ), and (3.4), one excludes amplitudes 
ipi and ipy , and finds that 



ip 2 ' = n {r) ^ 2 



where 



p* 



p 



(3.6) 



Eq. ( 1.12] ) together with Eqs. (3J3), and (|3.6| ), allows one 
to find the Green's function at any point of the loop. 

For instance, at point 2' and 2 where x^ r ' — and 
x (r) — respectively, 



1 



T pW 



(0 



R l-pU( l ) 



W 2 = P U^ 1 - - 



(0 



(3.7) 



(3.8) 



where T = 1 — i?. 

Although not obvious from the beginning, the struc- 
ture of the equations allows one to absorb the phase of 
p to the phase factors e lpF-C< , so that p becomes effec- 
tively a real positive parameter 



VR , R=\pf 



Then Eqs. (3/7), and (3J5) can be presented as 



U 2 = ^/RU {r) -TU {r} 
Uy = VRU (r) -T 



1-VflMC) ' 



i -Vrum 



U {r) 



Inserting the se ex pressions for the full evolution ma- 
trices into Eq. (2.6a), one gets after some transformation 
the Green's function g + at the points 2 and 2': 



9+2 = 1 + T 



.9+2' = 



I / Z/W 

1 + T 



1 -VflWM 



1 - 



l((r) 



These expressions for the Green's function are exact 
within the Andreev approximation where the normal 
metal spectrum is linearised in the vicinity of the Fermi 
surface; by a different method, similar formulae for a 
sandwich have been previously derived inU. In particu- 
lar, they allow for the space quantisation due to the nor- 
mal reflection (in addition to the AndreaMplevels grasped 
by the conventional quasiclassical theoryEra) . The evolu- 
tion matrices defined by Eq. Q, can be evaluated 
in the standard quasiclassical scheme, for instance, in the 
Riccati equation technique (see SectjBj). 



A. Checking the normalization condition 

Now, we want to demonstrate that unlike the simple 
loop case, the normalization condition for the Green's 
function for a double loop is violated even on the coarse- 
grained level. 

As before, the coarse grain averaging amounts to the 
integration of the Green's functions with respect to the 



(3.9) 
(3.10) 



i 

phases factors 



(3.11) 

£(i, r ) being the length of the corresponding loop. In prin- 
ciple, the integration-can be performed in terms of the 
elliptic integrals as inU. Instead, we use an approximate 
procedure which allows one to get better control of the 
physics behind the vio latio n of normalization. 

Analogously to E q. (|2.8 | ) , we present the evolution ma- 
trices entering Eq. (|3.9|) as 



u 



where u''' 1 *-* are given by Eq. (2.9) for each of the loops 



Algebra is greatly simplified if one uses the matrix rep- 
resentation where the evolution in the right (or left) loop 
is diagonal, i.e. 



y(r) _ e ip F c ir) 











(3.12) 
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This can be always achieved by the transformation of all 
the matrices, 



X -> oxer 



(3.13) 



with the matrix O, which diagonalizes U^ r \ In this rep- 
resentation, the right loop is effectively a normal metal. 
Up to a normalisation factor, the eigenfunction of 
can be written as: 



i (l) 



In this representation, A± has a simple physical meaning: 
A + (A-) are amplitudes of the Andreev reflection of the 
electron (hole) excitation for the ideal case t = 1, p = 
(a totally transparent interface). 

The left layer evolution matrix reads 

U <.»=e'.-'»(i(l +s ?>)«0) + i(l-rf>)^) 

(3.14) 



where 

„W - 



1 - A + A^ 



1 + A + A- 
2A+ 



-2A 
-l-A+A 



We plug these relations into Eq. 
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expansion in series of u^ 1 '^ . From Eqs. fl3.1^) , and (3.15) 
one sees that each power of u < - l ' r ' ) brings with it e +l(p 



(3.15) 
and perform an 



(electron) or e 



(hole). Accordingly, the expansion 



for a Green's function G has the following structure 



k n 



g= J2 (« (,) )*(« (r) ) n J2 

p=~k q——n 

(3.16) 



k,n=0 



Each term which can be associated with the correspond- 
ing process where the particle starts moving from an 
initial point and returns back having moved around the 
loops. En route, it undergoes electron-hole conversion as 
a result of the Andreev reflection on the knot. The Fermi 
momentum phase increases when the particle moves in 
the direction of the arrows in Fig. ^ (as the electron) 
and decreases otherwise (for the hole). The total ac- 
quired phase is controlled by the history of the electron- 
hole convers ion. The meaning of the integers k, n,p, and 
q in in Eq. ( 3.16 ) is the following. For instance, k shows 
the total length of the path kO- 1 ' in the left loop. If 
the particle have circled the loop k e times as electron 
and kh times as hole (fc = fc e + kh), the a cquir ed Fermi 
momentum phase corresponds to p in Eq. (3.16) equal to 
p = k e ~kh- Parameters n and q have analogous meaning 
for the right loop. 

On the coarse-grain level of description, one performs 
averaging within an interval of the loop length SC <C £ 
which is wide enough to average the phase factors in 



Eq. ( pi] ) to zero. Independent averaging with respect 
to length of the left and right loops, amounts to ignoring 
any term in Eq. (3.16) other than p = q = 0. For cer- 
tain geometries, the phase factors e l{pl ) and e IV< } cannot 
be considered independent. This occur 

if £W = £('') 

or 

more generally, the lengths of the loops are commensu- 
rable, iV£W = M£W , N,M = 1,2.... In this case, 
the terms p — ±N, q — =fM (and their multiples) also 
contribute to the coarse grain average. 

For dcfmitcncss, we consider the Green's function, <? 2 j 
corresponding to the point 2 in the immediate vicinity 
of the knot (see Fig. ||). The loops are assumed to be 
incommensurable. The result of the calculation of the 
averaged Green's function (32) can be presented as fol- 
lows: 



Gn 



(3.17) 



where G qc is of the quasiclassical Maki-Schopohl form, 



G qc — 



1 f 1 + a/3 -2a 

-a/3 ^2/3 -(l-a/3) 



with 



P = 



A +1 T ((l-^pvTu^ 2 

2 M 2 



((1 - 7 V^ ( ' )2 ) (1 + i 2 ^ 2 Ru {r)2 ) 

+ T 2 7 V )2 u «y (z/ 2 + ARA_ A+) ) + . (3.18) 



where 

/i = 1 - A + A- , v = 1 + RA + A- , 7 



1 - RA+A- 



these combinations of the Andreev amplitudes aws well- 
known in the theory of an imperfect NS-interfacolJ. The 
second term in the r.h.s. of Eq. (3.17), G nqci reads 



A+A- 
A + j{2TA + A- - pv) 





-A, A_ 



Squaring the coarse-grained Green's function, one gets 



(.92 



1 = 8A + A- p, 2 1 4 RT 2 u^ 



2 W 2 i 



(3.19) 



These expressions are valid up to u^ 2 and itW 2 ; their 
product included. 

One sees that the normalization condition is not re- 
stored by the coarse-grain averaging as in the cas e of a 
single loop. From the structure of the r.h.s. ofEq. ( 3.19 ), 
one may infer that the process responsible for the viola- 
tion of the normalization condition is the one where the 
particle follows an oo-like path twice along the loops, each 
in the opposite direction; such a path is (u^u^) 2 long 
and the change of the direction of motion requires two 
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sequential Andreev reflections, A + and A-. In the limits 
R — ► or T — ► 0, one restores the single loop situation 
and the normalization condition. Note also, the normal- 
ization is not violated if A± — 0, i.e. the electron- hole 
conversion is absent (in the normal metal case or if the 
loops are made of the same material) . 



IV. SANDWICH 

Now we are in position to apply the theory of previous 
sections to a physical system. We consider a two-layer 
system shown in Fig.[l] formed by two superconducting 
layer, left and right. The properties of the layers are 
specified by their thicknesses a('' r ) and the complex order 
parameter A"' r ) (they are considered as input, and the 
self-consistency condition is not discussed). We assume 
that the order parameter depends only on the coordinate 
x perpendicular to the layer plane. 

The Gor'kov equation for the matrix Green's function 
G e (ri, r 2 ) reads 



(e-H- 



= 1 



where HI is defined in Eq. (1.2) and V is the potential 
describing the the interface. As before, we assume the 
parabolic spectrum, i.e. 

"9 2 

| = P Pf 
2m 2m 

Due to the in-plane translational symmetry, the 
Green's function can be presented as 

G e (r,r') =J2 G ^ x > x ''>P\\y Pll ' (r ~ r ' )l1 ■ 
pii 

To simplify notations, we use below G( short- 
hand for G e (x,x';pu). 

The matrix 2x2 Green's function G(x, x') is found 
from 



The boundaries at x 
etrable so that 



are assumed impen- 



G(xi,X2 = 



= _„(') 



'/" or <i-' ') 



(ai = -a W or a {r \x 2 ) = 



(4.2) 



The Green's function G(x, x') is continuous at x = x' 

G(x + 0, x) = G(x - 0, x) , 

whereas its derivatives suffer a jump generated by the 
delta-function in the r.h.s. of Eq.(l.l): 



ill 
2m 



pGf z 



= 1 



If one chooses a particular model for the interface be- 
tween the layers, e.g. a (^-function barrier, one can solve 
the Gor'kov equation for the retarded Green's function 
exactly. The analysis of the full solution shows that the 
Andreev approximation, greatly simplifying the algebra, 
faithfully preserve all the physics we are going to discuss. 
Therefore, we restrict ourself to the Andreev approxima- 
tion as in Sect. | and III. 

As in Se ctio n ll we split the Green's function into two 
pieces Eq. (1.6), G + and G - , which describe propagation 
of the particle injected in the system with the momentum 
+PFx and —pFx, respectively. Consider for for definite- 
ness the Green's function G + in the right layer. Since the 
initial direction of the momentum may be changed by the 
reflection on the outer walls and the interface, instead of 



Eq. (1.8) one has 

iv F G + (xi,x 2 )t z 



g ++ {x u x 2 ) e <w.(*i-*3) 
9 + -{xi,x 2 ) e -^'« I1+l! - 2 " ,r ('l?) 



the function g ++ suffers the jump Eq. ( 1.1C| ) at the source 
point x\ — x 2 whereas the secondary wave is con- 
tinuous. 



From the boundary condition Eq. (4.2) 



,++, 



,0) 



,x 2 ) = g + {xi 



a ir) ,x 2 ) 



(4.4) 



where 



and 



(e-W x )G(x,x') = lS(x-x') 



e-i x -A 
-A* £ + 4 



Y(x) 



(4.1) 



( x = — (-4-X- P ^ , ^ = ^cos 2 , 
2m \ i dx I 2m 2m 2m 



where 



Because of the reflections, the particle moves along the 
trajectories shown in Fig. ^. 



(1) 



I 

-eJo— 



(r) 



— olo 

0,0 



FIG. 7: The trajectory coordinate. In each of the layers 
the trajectory coordinate x"' and are counted from the 
interface along the path of a particle which moves from the 
interface and is reflected on the surface boundary. 



cos 2 0=1 



Pl 



It is convenient to introduce the trajectory coordinates 
x( ' r >, which assume values in the region < x^ r > < 
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, as 



2a^ - x 
2a<' 



-fx 



< x r < a M ; 
< x < 2a( r ' . 

< x' < aW ; 
aW < x < 2a« . 



Since the trajectory coordinate space is twic e as large as 
the x— space, the two functions of x in Eq. ( |4.3| ) may be 
packed into a single function of the trajectory coordinates 
Xi and X2, that is 



IV F MXl,X2)T£ 



g + +(x 1: x 2 )e^^-^\ Xl ^ 2 
ff + -(a; 1 ,a ;2 )e-^((^-^- 2Q(r) )| Xl=2oM _ Xl!X2 



< xi < ; 
< xi < 2a( r ) ; 



for X2 < a,( r '. Similar relations, express &~ (x\, x 2 ) via 
( G(xi,x 2 ) for the value of x 2 in the region < x 2 < 
2a,( r \ The trajectory in the left loop is introduced in the 
same manner. 

Finally, one comes to the representation identical 
to t hat introduced for the loop Green's function by 
Eq. (O), i.e. 



:xi,x 2 )r 2 = 9 ( Xl ,x 2 )e^( xl - X3 ) 



IV F 

By virtue of Eq. (4.4), the quasiclassical envelope 
<?( x i! x 2) is continuous at the point of the loop corre- 
sponding to the outer walls. The envelope satisfies the 
quasiclassical equations Eq. (1.12) with the understand- 



ing x\^2 is the trajectory coordinate Xi. 2 . 

The coupling of the left and right layers can be de- 
scrib ed p henomenologically with the help of the S- matrix 
Eq. (3.1). The simple structure of the S-matrix in 
Eq. (3T) corresponds to the assumption the interface is 
mirror symmetrictS. In this context, p is the amplitude 
of the reflection by the interface and t is the transmission 
amplitude. Then, the theory of the loop Green's function 
presented in Sections |l| and III can be directly applied to 
the case of a layer or a sandwich. Within then Andreev 
approximation, this gives the solution which is "exact"jii 
the sense that it allows for the space quantization (incil, 
this Green's function is called quasiclassical-like). 

Schematically, the loops are shown in Fig. @: Apart 
from the Fermi momenta phase factor, the propagation 
of the waves moving towards the interface is described by 
the evolution operator and in the left and right 
layer, respectively. The evolution of the outgoing waves 

(I) (r) 

is controlled by Uy and f7 2 ' . The m " anc ^ out-waves 
are mixed by the S-matrix, an electron-hole scalar in the 
Andreev approximation. 



T he G reen's function can be then found from Eq. (3.E) 
and 3.1C where 

l/(l,r) — e (2pFxa <,,r) i+7ri) jj(l,r) 



(4.5) 



The evolution matrix depends on the direction of propa- 
gation in an anisotropic superconductor and for this rea- 



son and U^' may differ from each other. The ad- 
ditional 7r-phase shift is due to the reflection on the wall 



7-M 





FIG. 8: Waves with definite in-plane momentum in a sand- 
wich. The wave is totally reflected by the outer waves. The 
interface is partially transparent. To show the two-component 
electron-hole character of the wave, we use the double lines. 
The evolution operators U's mix electron-hole components 
whereas the interface scattering is described by an electron- 
hole scalar S-matrix. 



(the minus sign in the r.h.s. of Eq. (4.3)). The matrices 
can be found solving usual quasiclassical equations 
(see Section [§]). 

In the simplest case of isotropic homogeneous super- 
conductors, the evolution matrices are simply 



U U,r) 



i-K-\-2ippa^ l ' r ^ I cos 6 



5" 



{l,r) 
Sbulk 



1 



«('.»■) 

where ffouik stands for the bulk Green's function and 

u (l,r) = e . F »« 



£ is defined by Eq. (2.7) (in the ballistic case, £ = 
^ 2 ~|A| 2 ). 

T he tr ansformation to the convenient representation 
Eq. ( 3.13 ) where one of the layers (the right one) becomes 
effectively the normal metal can be performed with the 
matrixy 



O 



1 



''bulk 



i bulk 
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where o£Jk and a buU 



are the parameters of the Maki- 

,0) 



Schopohl parametrisation of S^uik 

It follows from the derivation in Sections || and [I], 
that the exact Green's function, e.g . at the poin t 2 (s ee 
Fig. g 2 can be found using Eq. ( ^6b| ) as Eq. ( |2~6b|) , 



ff2 



1 + 7L^U [r 
1 - TZ^U^ 



where in accordance with Eq. (3.6) TZ^ can be expressed 
via the reflection amplitude p and the evolution matrix 

As alr eady mentioned, the exact propagator given by 
Eq. ([[I]) must not be compared with the quasiclassical 
Green's function because they are different physical ob- 
jects. In deed , the Green's function in the Gor'kov equa- 
tion Eq. (4.1) describes a one-dimensional propagation of 
a plane wave with definite parallel momentum p|| and, 
therefore, of infinite extension in the direction parallel 
to the sandwich plane. Being infinite, the multiply re- 
flected waves inevitably overlap and their interference 
leads to the bound state formation (which would be re- 
veal itself as Fabry-Perot type resonances if the sandwich 
was connected to the outside world). On the other hand, 
a beam-like wave packet built of the plane waves with 
close values of pu, which behaves like a classical particle 
in a billiard, does not overlap with itself and the reso- 
nances are absent. The wave packet construction, which 
produces the trajectory picture shown in Fig. |[ is equiv- 
alent to the coarse-grained averaging discussed before. 
The coarse-gra ined Green's function for the sandwich is 
given by Eq. (3.17) at the point 2 chosen as a repre- 
sentative point. In what follows we want to show that 
in an inhomogenious superconductors the interference ef- 
fects survive the coarse grain averaging. For this, we find 
the quasiclassical Green's function and compare it with 
exact coarse-grained solution Eq. (3.17). 



where S is the interface S-matrix Eq. (3.1), and a = 
diag(ai,02) and b = diag(£>]/, b 2 i) are diagonal matrices 
bulit from the "input" Andreev amplitudes at the tra- 
jectory points immediately ajacent to the interface (Fig. 
g). The determinant reads 

V = R(l - aib v )(l - a 2 b 2 ,) + T(l - ai6 2 /)(l - a 2 b v ) , 

where R and T are the reflection and transmission coef- 
ficients. The "output" amplitudes ay^ and 61,2 can be 
found by the following formulae: 
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A. Quasiclassical Green's function 

In the framework of the quasiclassical theory, an it- 
eration procedure which allows one to find the Green's 
function of a sandwich, has been suggested inQ. Here, we 
use a modification of the method which is more conve- 
nient for the expansion with respect to the path length. 

We use the Maki-Schopohl parametrisation of the 
Green's functiorc3, 



1 



fJ 



1 



ab —2a 
1 - ab \ 2b -1 - ab 



and thc^bmindary conditions for the Andreev amplitudes 
a and ffl'trtl at the interface. 

Consider a collision with the interface, see Fig.||. In 
accordance withtJtil, the boundary condition can be con- 
veniently expressed via the determinant 

P = det||l-5'aS' t S|| 



FIG. 9: Scattering on the interface. The incoming channel are 
1 and 2, and the outgoing ones are 1' and 2'. The numbering 
is same as in Fig. H. The arrows show the direction of the 
momentum. The input for the interface boundary conditions 
is the set of Andreev amplitudes 01,2, &i',2'- These parameters 
can be found from the output of the scattering event on the 



previous collision a°i' 



■ old 
°1,2- 



In the iteration procedure illustrated by Fig. the 
updated input amplitudes, a\j, and by 2', are found from 



the output amplitudes of the previous step, 



,old 

n',2' 



and 



b° l 2, from the following relations: 



01,2 
1 



.old 



1 

h',2' 



where the "bond" evolution operators U^ 1 ^ are those in 
Eq. ([4.5|) . These relation are understood in the following 
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sense: For instance, to update a\ from the value a°] d , 

one calculat es th e column (°V ) , with the evolution 
matrix Eq. (3.14), and the advanced value a\ is found 
as the ratio of the upper and lower components of the 
column. 

Given the up date d input, the new output values are 
found from Eq. ( |4.6| ). The initial input is the bulk value 
of the parameters. In the following, we use the repre- 
sentation diagonalizing t he bu lk Green's function in the 
right layer (as in Section III A ) . Then the zero step value 
are b\ = A +1 ay = A-, 62 = wv = 0, a nd th e evolution 
operators are those in Eqs. ( 3.12| ), and ( 3.14L 

Unlike the iteration procedure suggested inl3, on any it- 
eration steps the parameters a and b have direct physical 
meaning - they give the Green's function at the central 
knot of the tree with 4" + 1 knots, n = 0, 1, . . . being the 
number of the iteration steps. 

The result of the calculation is the following. In the 
lowest order with respect to the path length, we recover 
the exact coarse-grained Gree n's fu nction discussed in 
Sect. |fv| and derived in Sect. Ill A. However, the qua- 
siclassical theory reproduces incorrectly terms which are 
proportional to the product u^ 2 u^ 2 , and, of course, 
higher order terms. Specifically, the quasiclassical the- 
ory faithfully re produ ces the first terms, G qc , in the exact 
expression Eq. ( |3.17 ) but it misses the second one, that 



is G nqc . We see that G nqc in Eq. ( 3.1 9| ) gives the 



cor- 



rection to the quasiclassical approximation in the lowest 
approximation with respect to the path length. 



reflected as a hole move b ack to "0' 
tion to u( r ) 2 -tcrms in Eq. (|3.18|)). 



(this give contribu- 
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FIG. 10: Typical trajectory in a sandwich, (a) The cross 
shows the initial position of the particle near the knot marked 
by "0" . Due to the possibility of the electron-hole conversion, 
the particle may change its branch and therefore the direction 
of motion when hitting the interface (dashed line). The con- 
nectivity of the trajectory is shown schematically in (b). The 
knots are numbered as in (a), letters r and I label the pieces 
of the trajectory in the right or left layer, respectively. This 
picture refers to the general case when the thickness of the 
left and right layers are incommensurable. Fig. ^(a) shows 
the case of a symmetric sandwich. 



V. DISCUSSION AND CONCLUSIONS 

In the previous section, we have shown explicitly that 
the quasiclassical theory gives, strictly speaking, incor- 
rect results when applied to a sandwich with a specular 
reflecting interface. By this, we confirm some previous 
observationstil. We understand this result as due to dif- 
ferences in the actual topology of the trajectories and 
those assumed in the quasiclassical theory. Our argu- 
ments are as follows 

In Fig.[l0|, we show fragment of the trajectory space 
which is large enough to include all the paths contribut- 
ing to the approximate Green's function in Eq. (3.17): 
the particle starts in the immediate vicinity of the knot 
"0" and after travelling around returns back. In the nor- 
mal metal case, the particle always follows the arrows 
being the electron, or always moves against the arrows 
if it is the hole, and there are no return paths. Corre- 
spondingly, coarse-grained properties are the same as in 
the bulk (as it is well-known from the the theory of nor- 
mal metal films). In a superconductor, the electron-hole 
conversion is possible due to the Andreev reflection pro- 
cess when the particle hits the interface separating the 
layers with different order parameters. Then, the parti- 
cle may move along the arrows as well as in the opposite 
direction. For instance (see Fig. [l0|(b)), the particle may 
move from "0" to "5" as electron, hit the "5" and being 



The connectivity in the quasiclassical theory is shown 
in Fig. O. The trajectory is tree-likeu where sequential 
collisions with the interface are completely uncorrelated 
as if the surfaces are rough: For instance, the knot 2% 
which can be reached from the knot 3, is regarded as 
different from the know 25 to which the particle may 
come from the knot 5. 
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FIG. 11: Tree like trajectory in the quasiclassical theory. 
Numbering of the notes in the same as in Fig. jj^. The 
subscript shows the knot from which the given knot can be 
reached. 

It is clear from Figs. and [II] that the whole class 
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of the paths allowed in the exact theory is not presented 
in the quasiclassical approach. These are loop-like paths, 
for instance, "0" "5" -> "2" -> "3" -> "0". This path 
includes two Andreev reflections - at the knots "5" and 
"3", and its length corresponds to u^ 2 vS r ^ 2 . The Fermi 
momentum phase factors cancel exactly (in a transla- 
tionally invariant system), and therefore, these paths are 
immune to the coarse-grain averaging. Summing over all 
the loop-like paths of this length, one is able to account 
for the r.h.s. of Eq. (3.19). In the geometry under con- 
sideration, the loops are not rare events (as in the weak 
localization theory), and, nominally, their contribution, 
i.e. the non-quasiclassical correction, is of order of 1. 

This analysis, which is the main result of the paper, 
gives the physical explanation of the invalidity of the qua- 
siclassical theory in many layer systems first noticed in 
Ref.Ej and analyzed numerically inO. The terms missing 
in the quasiclassical theory can be definitely attributed 
to the loops formed by the interfering paths. Similar 
conclusion has been made for the case of open geome- 
try calculating the Josephson current through a layered 
systemEa. The interference survives coarse-grain averag- 
ing because the lengths of the two intefering paths re- 
main equal for any initial condition. In turn, the loops, 
as closed orbits in the billiard theory, are possible only 
in very special - integrable - geometries. Here it is im- 
portant that the sandwich has a perfect in-plane homo- 
geneity and the reflection is specular. Any imperfection 
on the scale larger than the Fermi wave length (but, 
perhaps, negligible small compared with the coherence 
length), transforms the highly connective trajectory in 



Fig. [Tc](a) into the tree in Fig. |TT|. One may say that the 
quasiclassical theory is stable relative to small variations 
of geometry and is applicable in most of the cases. How- 
ever, there exist physically realizable situations where the 
conventional quasiclassical theory is inapplicable due to 
the contribution of the intefering paths. In this case, one 
should use the modification of the quasiclassical tech- 
nique suggested in the present paper. 

Finally, we want to mention that the physics due to the 
loops leads to a qualitative effect which is completely be- 
yond the quasiclassical theory. It is well-known that the 
gauge-invariant potential <I> does not enter the equation 
for the quasiclassical stationary retarded Green's func- 
tion, and therefore, the spectrum found from then Eilen- 
berger equation is completely unaware of its presence. 
However, the potential $ changes the actual value of the 
Fermi energy at which the Cooper pair condense and, 
therefore, modifies the value of the Fermi momentum by 
SpF = <&/vf- Varying one controls local value pf 
and the phase acquired by the particle on the intefer- 
ing paths. Proper choice of conditions, may turn the 
constructive interference into destructive, and, therefore 
modify the loop contribution to the density of states. 
Thus, the loop-like paths provide a physical mechanism 
for influencing the density of states in a superconductor 
by the non-equilibrium potential $. 
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APPENDIX A: OBSERVABLES 

For calculating observables the Green's function at co- 
inciding arguments G(x, x) is needed. It follows directly 
from the continuity of G(x, x') for x = x' and the defini- 
tion of g^' that 

G(x,x) = ^(d(+x,x) + G(-x,x)j = 

= i (G (+) (+x, x) + G { - ] (+x, x) + G (+) (-x, x) 
= \ (g W {x)+g(-\xj) . 



function g(x) 
1 



2 tanh— &Tr 



(g {+ \x) - g(-\x))r z 



The local density of states can be written in terms of G 
as follows 



= — lim Tr (G(x, x') - G A (x, x')) . 

Z7T x— *x' 

From here the following expression for LDOS via the 1- 
point Green's function g^ (x) can be found 



v(x) = 



2ttvf 



-KTr 



g^(x)+g(-\x) 



This is the formula for the local density of states at a 
point x expressed in terms of the 1-point Green's function 

The standard quantum mechanical expression for the 
current density via Gor'kov Green's function reads 

In the equilibrium the Keldysh Green's function can be 
expressed as <G K = (G- G A )tanhe/2T. The differenti- 
ation with respect to x and x' is performed using identi- 
ties d x G& = ±i Pa d( ± ) and d x ,G^ = TipcG^. This 
leads to 



-v F 



tanh -J- lim 3? — Tr f G (+) - G^ A 

2T x— >x' Wf \ / 



Finally we take the limits x — > x' and get the expres- 
sion for the current density in terms of 1-point Green's 



APPENDIX B: EVOLUTION MATRIX 



Q(-3?he j)vn-pV)se of this section is to present a method for 
calculatlonAhe quasiclassical evolution matrix. 

One can always present U(x,x') in the following form 

U(x,x') = a(x,x')( . 1 ^ 0(1,0) + b(x,x')(^ X ' ^ <g> (0,1) 
\a{x,x')) \ 1 J 



where 



(Bl) 



a(x, x) = b(x, x) = 1 , a(x, x) = (3(x, x) = (B2) 

as required by the condition U(x, x) = 1. Since the evolu- 
tion matrix must generate a solution fo r ar bitrary initial 
condition at x', the two terms in Eq. (Bl) must satisfy 
the quasiclassical equation separately. 

From Eq. ( |2.3| ). the two-component "wave function" 
(") = ( a ) an( i = ( h ) considered as functions of x, sat- 



isfy the Andreev-type linear first-order.d inferential equa- 
tion. It is known from the literatureEHtJ that the ratio 
u/v obeys the nonlinear Riccati equation, 

i— a = 2ea + A* R + Aa 2 , -i— f3 = 2e/3 + A* R f3 2 + A 
ox ox 



with the initial condition in Eq. (|B2|). 

Known o_and /3, parameters a{x,x') and b(x, x 1 ) can 
be found astl 



i J dx(e+Aa) -i f dx(s+A* R p) 

a(x,x') = e*>' , b(x,x') = e x ' 



This completes the derivation of the evolution matrix in 
a general case. 



